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Abstract 



Aims. The comparative study of several molecular species at the origin of the gas phase chemistry in the ditfuse interstellar medium 
' ' (ISM) is a key input in unraveling the coupled chemical and dynamical evolution of the ISM. 

' Methods. The lowest rotational lines of HCO+, HCN, HNC, and CN were observed at the IRAM-30m telescope in absorption against 

the A3 mm and Al.3 mm continuum emission of massive star-forming regions in the Galactic plane. The absorption lines probe the gas 
. over kiloparsecs along these lines of sight. The excitation temperatures of HCO^ are inferred from the comparison of the absorptions 

• in the two lowest transitions. The spectra of all molecular species on the same line of sight are decomposed into Gaussian velocity 

components. Most appear in all the spectra of a given line of sight. For each component, we derived the central opacity, the veloc- 
ity dispersion, and computed the molecular column density. We compared our results to the predictions of UV-dominated chemical 
?^ , models of photodissociation regions (PDR models) and to those of non-equilibrium models in which the chemistry is driven by the 

dissipation of turbulent energy (TDR models). 
' Results. The molecular column densities of all the velocity components span up to two orders of magnitude. Those of CN, HCN, and 

HNC are linearly correlated with each other with mean ratios A'(HCN)/A'(HNC) = 4.8 ± 1.3 and A'(CN)/iV(HNC) = 34 ± 12, 
and more loosely correlated with those of HCO+, MHNO/MHCO^) = 0.5 ± 0.3 , A'(HCN)/A'(HCO+) = 1.9 ± 0.9, and 
A'(CN)/A'(HCO^) = 18 ± 9. These ratios are similar to those inferred from observations of high Galactic latitude lines of sight, 
^ ' suggesting that the gas sampled by absorption lines in the Galactic plane has the same chemical properties as that in the Solar neigh- 

CNj I bourhood. The FWHM of the Gaussian velocity components span the range 0.3 to 3 km s"' and those of the HCO^ lines are found to 

00 . be 30% broader than those of CN-bearing molecules. The PDR models fail to reproduce simultaneously the observed abundances of 

l/^ ■ the CN-bearing species and HCO^, even for high-density material (100 cm"^^ < < 10"* cm"''). The TDR models, in turn, are able 

' to reproduce the observed abundances and abundance ratios of all the analysed molecules for the moderate gas densities (30 cm"^ 

< «H < 200 cm"^) and the turbulent energy observed in the diffuse interstellar medium. 

Conclusions. Intermittent turbulent dissipation appears to be a promising driver of the gas phase chemistry of the diffuse and translu- 
cent gas throughout the Galaxy. The details of the dissipation mechanisms still need to be investigated. 
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1. Introduction tre, and centimetre wavelengths in absorption against the contin- 

uum emission of star-forming regions (e.g. Koo 1997, Fish et 
.C^ Smce Its discovery through absorption lines m bright star spec- 2003, Nyman 1983, Nyman & Millar 1989, Cox et al. 1988, 
tra (Hartmann 1904), our knowledge of the interstellar medium ^^^^^ ^ ^^j^j^ ^992, Greaves & Williams 1994, Neufeld et 
has grown thanks to a variety of absorption hne measurements. 2002, Plume et al. 2004, Olofsson et al. 2010 in prep.) and 
Diatomic molecules have been discovered in this hostile en- ^^- ^^^ extragalactic radio sources (Liszt et al. 2008 and refer- 
vironment in the late 1930 s and eai-ly 1940 s (see references ^^^^^ ^^^^.^-^^ -^^^^.^ ^^^^ ^^^^ ^^^^^ measure- 
in the review of Snow and McCall 2006) and since then have ^^^^^ complex, and the link between the structure in density 
been observed in the UV and visible specti-al domains toward temperature, the molecular richness, and the velocity field 
bright stars (e.g. Crane et al. 1995, Gredel 1997, Weselak et ^^^^ ^^^^^ unraveled. In particular, the abundance of several 
al. 2008a, 2008b, 2009, Gry et al. 2002, Lacour et al. 2005) molecules, such as CH+ or HCO+, are found to be at least one or- 
at increasingly high extinction values (up to Ay ^ 5, Gredel et magnitude larger than those deduced from UV-dominated 
al. 2002). Molecules are a lso detected at submillimetre, millime- chemical models. In the cases of HCO+, HNC, and HCN the 

J, _ ' ; '. , , ■ , . Tr, » , J 1 abundances inferred from observations of the diffuse ISM are 

Based on observations obtained with the IRAM 30m telescope. ... , , n r • 

IRAM is supported by INSU/CNRS (France), MPG (Germany), and similar to those observed in dark clouds (Lucas & Liszt 1996, 

IGN (Spain). 2000, Liszt & Lucas 2001). The column densities of specific 
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Table 1. Properties of background sources, and rms noise levels of the spectra. 



Source 


RA(J2000) 


Dec (J2000) 


1 


b 


D" 






o-,IT, " 


a-, IT, " 


0-,/r, " 




(h) (m) (s) 


n (') (") 


O 


n 


(kpc) 


HCO+ (0-1) 


HCO+ (1-2) 


HNC (0-1) 


HCN (0-1) 


CN(O-l) 


G05. 88-0.39 


18 00 30.4 


-24 04 00 


05.88 


-0.39 


3.8 


0.021 










G08.67-0.36 


18 06 18.9 


-21 37 35 


08.67 


-0.36 


4.8 


0.138 










GlO.62-0.38 


18 10 28.7 


-19 55 50 


10.62 


-0.38 


4.8 


0.049 


0.206 


0.043 


0.035 


0.064 


G34.3-I-0.1 


18 53 18.7 


+01 14 58 


34.26 


+0.15 


3.8 


0.117 


0.328 


0.038 


0.036 


0.154 


W49N 


19 10 13.2 


+09 06 12 


43.17 


+0.01 


11.5 


0.015 


0.072 


0.020 


0.021 


0.031 


W51 


19 23 43.9 


+ 14 30 30.5 


49.49 


-0.39 


7.0 


0.022 


0.089 


0.043 


0.034 


0.070 



" Source distance from Fish et al. (2003) wiio resolve the kinematic distance ambiguity. Errors are around 1.0 kpc. 
* (Tz/Tc is the rms noise divided by the continuum intensity of the spectra, at a frequency resolution of 40 kHz. 



Table 2. Observation parameters. 



Molecule Transition vo (GHz) (K) T^ys (K) 

G05 G08 GIO G34 W49 W51 G05 G08 GIO G34 W49 W51 



HCO+ 


0- 


1 


89.1885247 


1.5 0.22 0.69 


0.86 


1 


73 


1.03 


156 162 206 


137 


167 


132 


HCO+ 


1 - 


2 


178.3750563 


0.60 


0.26 


1 


10 


0.96 


588 


420 


375 


374 


HNC 


0,1 - 


1,0 


90.6634170 


0.97 


0.88 


1 


94 


0.93 


186 


162 


160 


191 


HNC 


0,1 - 


1,2 


90.6635560 


0.97 


0.88 


1 


94 


0.93 


186 


162 


160 


191 


HNC 


0,1 - 


1,1 


90.6636220 


0.97 


0.88 


1 


94 


0.93 


186 


162 


160 


191 


HCN 


0,1 - 


1,1 


88.6304160 


1.05 


1.20 


1 


86 


1.10 


158 


167 


146 


140 


HCN 


0,1 - 


1,2 


88.6318470 


1.05 


1.20 


1 


86 


1.10 


158 


167 


146 


140 


HCN 


0,1 - 


1,0 


88.6339360 


1.05 


1.20 


1 


86 


1.10 


158 


167 


146 


140 


CN 


0,1/2,3/2- 


1,1/2,1/2 


113.1441573 


0.94 


0.70 


1 


68 


0.92 


280 


382 


275 


275 


CN 


0,1/2,1/2- 


1,1/2,3/2 


113.1704915 


0.94 


0.70 


1 


68 


0.92 


280 


382 


275 


275 


CN 


0,1/2,3/2- 


1,1/2,3/2 


113.1912787 


0.94 


0.70 


1 


68 


0.92 


280 


382 


275 


275 


CN 


0,1/2,1/2- 


1,3/2,3/2 


113.4881202 


0.93 


0.69 


1 


68 


0.93 


284 


312 


209 


206 


CN 


0,1/2,3/2- 


1,3/2,5/2 


113.4909702 


0.93 


0.69 


1 


68 


0.93 


284 


312 


209 


206 


CN 


0,1/2,1/2- 


1,3/2,1/2 


113.4996443 


0.93 


0.69 


1 


68 


0.93 


284 


312 


209 


206 


CN 


0,1/2,3/2- 


1,3/2,3/2 


113.5089074 


0.93 


0.69 


1 


68 


0.93 


284 


312 


209 


206 


CN 


0,1/2,3/2- 


1,3/2,1/2 


113.5204315 


0.93 


0.69 


1 


68 


0.93 


284 


312 


209 


206 



species, like HCO^ and OH, exhibit remarkable correlations that 
cannot be understood in the UV-dominated chemistry. 

Alternative chemical models that couple the chemical evo- 
lution of the gas to the turbulent dynamical evolution of the 
medium have been developed. The space time intermittency of 
turbulent dissipation is invoked to locally enhance the rate of 
highly endoenergetic reactions, otherwise blocked in the cold 
ISM. Turbulent dissipation in low-velocity, magneto hydrody- 
namical (MHD) shocks (Flower & Pineau des Forets 1998) or 
magnetized regions of intense velocity shear (Falgarone et al. 
1995, Joulain et al. 1998) are promising frameworks. Turbulent 
mixing between the warm neutral medium (WNM) and the cold 
neutral medium (CNM) has also been proposed as a possibil- 
ity to enhance the formation rate of specific species (Lesaffre et 
al. 2007). In the turbulent dissipation regions (TDR) model of 
Godard et al. (2009), dissipation of turbulent energy occurs in 
short-lived (~ 10^ yr) magnetized vortices and is responsible for 
long-lasting (~ 10^ yr or more) chemical signatures due to the 
strong thermal and chemical inertia of the diffuse gas. Following 
the chemical and thermal evolutions of the dissipation and relax- 
ation phases, the TDR models reproduce the column densities of 
CH+, CH, HCO+, OH, H2O, C2H, and of the rotational levels of 
H2 (J > 3), as well as their correlations, observed in the local 
diffuse medium. 

The present study broadens the investigation of the local dif- 
fuse medium chemistry by analysing observations performed 
in the direction of remote star-forming regions that sample 
gas throughout the Galactic plane. The selected background 



sources are extensively studied star-forming regions, close to the 
Galactic plane {\h\ < 0.5°) (e.g. Mookerjea et al. 2007). These 
lines of sight are also the primary targets of the Herschel-HIFI 
(Heterodyne Instrument for the Far Infrared) key programme 
PRISMAS (PRobing InterStellar Molecules with Absorption 
line Studies) whose objective is to advance the understanding 
of astrochemistry by observing species at the origin of gas phase 
chemistry, namely light hydrides and small molecules containing 
carbon. The present work is therefore anticipating future com- 
parisons between radio and far infrared (FIR) observations. 

The observations and their analysis are discussed in Sects. |2] 
and|3] The main results are presented in Sect. |4] and the cyanide 
chemistry and the predictions of ID chemical models are dis- 
cussed in Sect. |5] 

2. Observations 

The observations were carried out at the IRAM-30m telescope 
at Pico Veleta (Spain) in August and December 2006. For the 
sources listed in Table [T] (with their Galactic coordinates and 
their distance from the Sun), we observed in wobbler-switching 
mode: 

• the y = - 1 absorption lines of HCO+ and HNC, 

• the F=l-l,l-2 and 1-0 hyperfine components of the 
7 = 0-1 absorption line of HCN, and 

• the F = 3/2 - 1/2, 1/2 - 3/2, 3/2 - 3/2, 3/2 - 5/2, and 
1 /2 - 1/2 hyperfine components of the 7 = 0-1 absorption 
line ofCN. 
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G1 0.62-0.39 G34.3+0.1 W49N W51 




-80 -40 40 80 -40 40 80 120 -80 -40 40 80 -80 -40 40 80 

LSR Velocity (km.s"^) 




10 20 30 40 50 10 20 30 40 50 30 40 50 60 70 -10 10 20 30 



LSR Velocity (km.s"^) 

Figure 1. Absorption profiles observed in the direction of GlO.62-0.39, G34.3H-0.1, W49N, and W51 in the ground state transitions 
of HCO^, HCN, HNC, and CN. The hyperfine structure (relative positions and relative LTE line strength) of the 7 = 0-1 transition 
of HNC and HCN are displayed. The broad velocity coverage of the upper panels illustrates the quality of the baseline and the 
complexity of the emission and absorption mixture. The lower panels display the velocity structure of the absorption features in 
more details. All the spectra of the lower panels have been normalized to the continuum temperature. 



We used the IRAM-30m SIS receivers tuned in SSB (single side quency resolution of 40 kHz. The rms pointing accuracy of the 
band) mode, with rejection of the image bands larger than 15dB. telescope checked on nearby continuum sources was 3" and the 
The spectra were obtained using the VESPA correlator at a fre- integration time ranged between 1 1 and 59 minutes. Between 88 
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and 113 GHz, the forward and main beam efficiencies were of 
0.95 and 0.75, the spectral resolution and the half power beam 
width range from ~ 0.13 km s"' and 28" to ~ 0.10 km s"' and 
22" respectively. The observation parameters are listed in Table 

12 

The gain receiver stability, and therefore the uncertainty on 
the continuum level, T^, is estimated by comparing the contin- 
uum levels from several spectra in the vicinity of the different 
lines. For the same frequency ranges, they differ by 0.1 to 0.3 K 
(see Table |2]l- The uncertainty on the continuum level therefore 
ranges between 10 and 30 %. 

Fig. [1] displays a selection of spectra obtained after data re- 
duction using the GILDAS-CLASS90 softwarf] (Hily-Blant et 
al. 2005). We focus on the absorption part of the spectra, outside 
the line emission of the star-forming region, since we are inter- 
ested in the velocity structure and properties of the absorbing 
gas. In several cases, in order to extract as much information on 
the absorption lines as possible, the emission in wings of strong 
lines from the background source was removed using polyno- 
mial and exponential fitting routines. 

As shown in Fig [1] the absorption spectra are highly struc- 
tured. While all the background sources are in the Galactic 
plane along lines of sight that cross the molecular ring and the 
Sagittarius spiral arm once (except that toward W49N which 
crosses it twice), the number of absorption features is highly 
variable from one source to another: ranging from 1 to 20 com- 
ponents. Optical depths, defined as 



are also highly variable from one species to another; several 
components are only detectable in HCO^, and the CN spec- 
tra are unequivocally less overcast than all the other spectra. 
Complementary data have been obtained with the PdBI interfer- 
ometer in the direction of W49N and W5 1 (Pety et al., in prepa- 
ration). 

The hyperfine components of the 7 = 0-1 absorption lines 
of HNC (Bechtel et al. 2006) are too close to be individually re- 
solved given the significant velocity dispersion of the gas. The 
hyperfine structure causes a systematic broadening that depends 
on the FWHM (full width at half maximum) Avi-eai of the ve- 
locity component: the broadening Ay - Awreai ranged between 
0.19 and 0.035 km s ' for Une FWHM varying from 0.3 to 3 km 
s~^ In the case of HCN and CN (0-1), the separation of the hy- 
perfine components is larger than most linewidths but the large 
number of velocity features along the lines of sight GlO.62-0.38 
and W49N induces a blending of the hyperfine components. 

While absorption features above the noise level appear at all 
widths, down to the velocity resolution, we have not analysed 
the spectra to that level: every feature we detect spans at least 3 
spectrometer channels, between 0.3 and 0.4 km s"' depending 
on the frequency. 

3. Analysis 

3.1. Excitation temperatures 

We observed the HCO^ (1-2) line toward 4 of the 6 sources 
listed in Table [T] As shown in Fig.|2 which displays the opaci- 

' See http://www.iram.fr/IRAMFR/GILDAS for more information 
about GILDAS softwares. 

^ This result on the line profile broadening is derived from the anal- 
ysis of 560 synthetic spectra taking into account the hyperfine structure 
of HNC (line strength and velocity structure). 



ties thco+(0-i) and thco+(i-2) in the absorption Unes channel by 
channel, the rms noise levels are low enough to compute ac- 
curate values of the excitation temperature toward W49N and 
W51. Toward GlO.62-0.39 and G34.3-H0.1, the correlation be- 
tween thco+(0-i) ™d Thco+(i-2) is looser and could be due to 
actual variations of the excitation temperature along the line of 
sight. 

However, for the sake of simplicity, we assume that the ab- 
sorbing gas on each line of sight is defined by a single excitation 
temperature Tex, derived as 

l-exp(-/iV2/fcrex) ^ J-ridv lv2? g,,iA»n 

eXp(/!Vi/ferex) - 1 jTidv\Vlj gulAull 

where ri and T2 are the line opacities per unit velocity, v\ and 
V2 the rest frequencies, gu\ and g„2 the upper level degenera- 
cies, and A„;i and A„/2 the Einstein's spontaneous emission co- 
efficients of the molecular transitions HCO^ (0-1) and (1-2) re- 
spectively. The results are listed in Table [3] For all sources, Tex 
is found close to the temperature of the cosmic microwave back- 
ground TcMB = 2.73 K. 

Table 3. Excitation temperatures inferred from the 

'!"hco+(i-2)/'''hco+(0-i) ratio. 



Source 


Tex (K) 


errors" 


(K) 






+ 




GlO.62-0.38 


3.0 


1.4 


0.5 


G34.3+0.1 


3.0 


8.4 


0.8 


W49N 


2.7 


0.8 


0.4 


W51 


2.9 


1.1 


0.5 



" The errors are derived from the results of numerical calculation 
which do not take into account the lower limit of 2.73 K set by the 
cosmic microwave background. 



These low excitation temperatures suggest that HCO^, HNC, 
and HCN, whose dipole moments are similar, are radiatively 
rather than collisionally excited. Using the Large Velocity 
Gradients (LVG) code by Schilke (private communication), we 
obtain an upper limit on the gas density: hh < 10'* cm""*. This 
constraint is not stringent because of the large critical densities 
of the transitions (~ 10^ and ~ 10'' cm"^^ for HCO+ (1-0) and 
HNC and HCN (1-0) respectively). 

Estimates of the gas densities are provided by the LVG anal- 
ysis of emission lines observed at the velocities of the absorp- 
tion components. Toward W49N, '-CO and '^CO(l-O) and (2-1) 
line observations provide H2 densities all close to 5 x 10^ cm"^ 
(Vastel et al. 2000). Including the CI ^Pi -^Pq Hne in the analysis 
of the CO lines. Plume et al. (2004) find lower densities, rang- 
ing between 1500 and 3000 cm"^. In the following, we therefore 
adopt hh < 5 X 10^ cm"^ as an upper limit of the gas density 
causing the absorption features. 

3.2. Decomposition of tiie spectra into Gaussian 
components 

The decomposition of the spectra in velocity components and the 
resulting column densities are inferred from a multi-Gaussian 
fitting procedure based on the Levenberg-Marquardt algorithm 
and the following sequence. 
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One of the absorption spectra is decomposed without a pri- 
ori into the minimal number of Gaussians required to fit the 
data within the observational errors. We preferentially use 
the CN (0-1) transition because its hyperfine structure pro- 
vides a valuable constraint on the line centroids. 
Based on the results of the previous step, the algorithm is 
then applied recursively to the others absorption spectra, rul- 
ing out a shift by more than a resolution element between 
the velocity components of each molecular species. In a few 
cases, a shift of at most 4 (~ 0.5 km s"') resolution elements 
is allowed because it ensures the convergence of the fitting 
routine: for example in regions of the spectra where many 
components are blended, as in the HCN (0-1) data at 59.5 
km s ' toward W49N. 



Thus, for each transition, 
(line/continuum) is written: 



the observed line profile 



' mb 



(v) = exp 



N, Ni, 

■ ^ ^ aiAk) To( j) e 

j=\ k=\ 



r ..-..Q(j)-i.o,,(t) i 



(3) 



where and are the main beam and the continuum temper- 
atures, A^e and Nh are the numbers of components and hyperfine 
transition^ tq, vq and Uv are the Gaussian parameters, and a/, 
and uo/i are the LTE and optically thin line strengths relative to 
a reference hyperfine transition and the velocity shift associated 
to each hyperfine transition. As described in Appendix lAl the 
column density associated to each velocity component is then 
derived assuming a single excitation temperature Tex for all the 
levels of a given molecule. 

For the components which are either saturated or too weak 
to be singled out by the fitting procedure, limits on the column 
densities are derived from the integrated optical depth over the 
corresponding velocity range. All the results are listed in Tabs. 
IA.HIA.4l of Appendix lAl and, as an example, the outcome of the 
multi-Gaussian decomposition, applied to the absorption lines 
observed toward W49N, is shown in Fig. ??. We find that the 
widths of the Gaussian components have a continuous distribu- 
tion with values ranging from 0.3 km s"' to 3.5 km s"', inde- 
pendently of the source or the molecular species. The peak opti- 
cal depths range between 0.06 and 2.2 and the inferred column 
densities per velocity component span more than one order of 
magnitude on each line of sight. 

3.3. Systematic errors on calculation of column densities 



The uncertainties given in Tabs. lA.lllA.4l are the formal 1-cr er- 
rors derived from the diagonal elements of the covariance matrix 
and do not take into account the systematic errors introduced by: 
(1) the finite velocity resolution, (2) the uncertainty on the exci- 
tation temperatures (see Table [3]), and (3) the uncertainty on the 
continuum level. 

We show in Appendix |B] that the finite velocity resolution 
introduces an error on the column densities smaller than 12 %. 
Table [3] gives the uncertainties on the excitation temperature 
of the lowest rotational levels of HCO^, that affects the parti- 
tion function, hence the column densities. Toward GlO.62-0.38, 
W49N, and W51, the uncertainties on Tgx are small, providing 
uncertainties on the column densities <^?aJ. Toward G34. 3+0.1, 
the inferred column densities could be underestimated by a fac- 
tor of 3. 



Finally, the uncertainty e = ST^/Tc on the continuum tem- 
perature (see Sect. 2) introduces an error 6t on the calculation of 
the optical depths: 



6t ~ e - \n(l + eeO . 



(4) 



This error is larger, in most cases, than those inferred from the 
fitting procedure, and for e = 10% the corresponding uncertainty 
on the column densities ranges from 16% to 62% when tq varies 
between 0.5 and 2. All components with tq > 2.2 are thus con- 
sidered as saturated. The components at 16 and 41 km s"' of the 
GlO.62-0.38 HCO+(0-1) spectrum are the only ones that fall in 
this category. 

When all uncertainties are taken into account, the column 
densities are determined within a factor of 2, and within a factor 
of 3 for the absorption lines observed toward G34.3+0. 1 . 

4. Results 

4.1. Comparison of column densities and profile imewidths 
among species 

Because the line centroids of a Gaussian component observed 
in several transitions coincide within 0.15 km s"' (except for 
6 components, see Appendix lAl Tabs. lA. ll - IA.4l l. we propose 
that the corresponding velocity component has a physical reality, 
in support of the comparison between the optical depths in the 
different molecular lines and the resulting column densities. 

The column densities and linewidths of all the unsaturated 
velocity components are displayed in Figs. |4] and ?? for sev- 
eral pairs of species. These figures also include, for comparison, 
the results of previous studies of molecular absorption lines ob- 
served toward strong extragalactic mm-wave continuum sources 
(Lucas & Liszt 1996, 2000 ; Liszt & Lucas 2001) that mosdy 
sample diffuse gas in the Solar Neighbourhood. 

Although we have not made any assumption on the nature 
of a Gaussian velocity component, it is remarkable that (/) the 
range of the component linewidths (0.3 to 3.5 km s"') is about 
the same in both data sets, and ( it) the range of column densi- 
ties of most molecular species (up to two orders of magnitude) 
is also similar in both sets. This result suggests that the gas com- 
ponents sampled at high galactic latitude and toward the inner 
Galaxy share common kinematic properties. In addition the rota- 
tion curve in the first quadrant provides displacements ~ 70-100 
pc/km s"' : therefore, if we assume that the induced spreading of 
the velocity centroids is responsible for the observed linewidths, 
a component of 3 km s"' (resp. 0.3 km s"') would have a size of 
300 pc (resp. 30 pc). This would correspond to densities «// < 40 
cm"^, below our estimates (see Sect. 14.3b . Hence, the linewidths 
of the components are likely due to turbulence rather than differ- 
ential rotation of the galactic plane. In the following, we com- 
pare the chemical and kinematic properties of these two ensem- 
bles. While the thick lines in Figs. |4] and ?? result from a linear 
regression unweighted by the fractional errors on the fitting pa- 
rameters, the values given below come from weig htecfl averages. 

• Concerning the column densities of HNC, HCN, and 
CN, Liszt & Lucas (2001) found tight linear correla- 
tions along lines of sight at high Galactic latitude, namely 
A^(HCN)/A^(HNC) = 4.8 + 1.1 and A?(CN)/A?(HNC) = 



A';, depends on the molecule 



A weighted average makes more sense numerically but tends to 
favour the high column density points. In practice, both evaluations 
(weighted and unweighted) are meaningful since the errors are not 
known with high accuracy (see Sect. l3.3t . 



6 



B. Godard et al: Molecular absoqjtion lines toward star-forming regions 



1e+13 




o 

^1e+12 



1e+11 



—I I I 



1e+10 1e+11 1e+12 1e+13 
N(HNC) (cm"^) 



0.0 1.0 2.0 3.0 
Av(HNC) (km.s"^) 



Figure 4. Comparisons of the column densities (left panels) and the FWHM (rigth panels) of HCO^, CN, HCN, and HNC. The red 
filled symbols are from this work. The black open symbols are from Lucas & Liszt (1996) and Liszt & Lucas (2001). The solid lines 
result from a linear regression of the data of the present work unweighted by the fractional errors on and Av : Y - a X. 



33 + 12. In comparison the mean ratios of the present work 
are A?(HCN)/A?(HNC) = 4.8 ± 1.4 and A^(CN)/A?(HNC) = 
34± 12. 

• Fig- S (top) shows that the coiTelations of the CN-bearing 
molecules with HCO^ are looser and probably non lin- 
ear with mean ratios: A^(HNC)/A^(HCO+) = 0.5 ± 0.3 , 
Ar(HCN)/A^(HCO+) = 1.9 + 0.9, and A^(CN)/A?(HCO+) = 
18 + 9. 

• Last, the loose linear relation between the column densities 
of C2H and HCO^ observed in Fig. ?? corresponds to a mean 
ratio A?(C2H)/A?(HCO+) = 8.3 ± 5.6 to be compared with 



A^(C2H)/A?(HCO+) = 14.5 + 6.7 obtained by Liszt & Lucas 
(2001). 

The linewidths of the Gaussian components of the HCO^, 
CN, HCN, and HNC line profiles are close to be the same. We 
find that the CN, HCN, and HNC line profiles are systematically 
narrower than those of HCO^ (by a factor of 0.7 to 0.9). In Fig.|4] 
the CN and HCN line profiles appear to be narrower than those 
of HNC (by a factor of 0.8 to 0.9) but this result is mostly due 
to the hyperfine structure of HNC that has not been taken into 
account in the decomposition. Once the coiTection discussed in 
Sect.|2]is applied, the significance of the linewidth difference dis- 
appears. The profiles of C2H and C3H2 have been studied with a 
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different method, and the components identified are not the same 
as those in the present study (Gerin et al. 2010). 

4.2. Estimation of fhe total liydrogen column densities along 
the lines of sight 

In order to be able to compare the data to chemical models, 
molecular abundances must be derived and column densities of 
hydrogen measured. In particular, it is essential to determine 
whether the wide dynamic ranges over which the correlations 
are observed (Figs.|4]and ??) are related to variations of the total 
column density of the gas sampled or/and of the physical and 
chemical conditions in the absorbing gas. The main difficulty is 
to estimate the fraction of molecular hydrogen that is not directly 
observable. Using All cm observations of HI, A9 cm and AQ.3 
mm observations of CH, and the remarkable correlation between 
CH and H2, A^(CH)/A?(H2) = 4.3 x lO"'* (Liszt & Lucas 2002), 
we evaluate the total amount of gas A^h - N(H) + 2N{H2) along 
the galactic lines of sight, as in Godard et al. (2009) for the lines 
of sight studied by Liszt & Lucas (2001). The HI column densi- 
ties are inferred from VLA All cm absorption line observations 
(Koo et al. 1997; Fish et al. 2003). Wherever possible, we derive 
A^(H2) from CH observations (at 9 cm by Rydbeck et al. 1976, 
at 1 THz by Gerin et al. 2010 in prep.). 

An independent estimate of the total column density of gas 
toward the star-forming regions is inferred from the analysis of 
the 2MASS survey (Cutie et al. 2003). Marshall et al. (2006) 
have measured the near infrared colour excess in large areas of 
the inner Galaxy (|/| < 100°, \b\ < 10°) to obtain the visible 
extinctions (Ay ~ lOA^:), providing an estimate of the total hy- 
drogen column density along the lines of sight. 

Table ID lists the HI (and H2, where available) column den- 
sities in selected velocity intervals, as well as the total hydro- 
gen column densities inferred from extinction. However the un- 
certainties on these estimations are large: (1) the error on the 
A^(CH) - A^(H2) relation is about a factor of 3 (Liszt & Lucas 
2002) ; (2) this correlation has been established in the local dif- 
fuse medium but has never been observed in the inner Galaxy 
material ; (3) because of the low resolution of the 2MASS survey 
(~ 15 arcmin), the error on the total hydrogen column density 
(computed as the standard deviation of the extinction measured 
along the four closest lines of sight surrounding a given source) 
is larger than 30 % ; and (4) The HI column densities inferred 
from VLA All cm absorption line observations are directly pro- 
portional to the assumed spin temperature. Hence, while the two 
determinations agree with each other within 15% toward W51 
and W49N, they differ by at least a factor of 2 toward G 10.62- 
0.38. 

According to the extinction measurements, the lines of sight 
sample between 1.3 (W51) and 12.5 (W49N) magnitudes of gas. 
The total velocity coverage of the absorption features is ~ 10 
km s~' toward W51 and 48 km s"' toward W49N. Therefore, 
the average hydrogen column density per velocity unit is only 
twice larger along the line of sight toward W49N than along 
that toward W5 1 . It is therefore possible to estimate the molecu- 
lar abundances relative to the total hydrogen column density A^h 
for each velocity component, assuming that A^h scales with their 
linewidth according to Nh/^v - 1.1 and 4.7x10^" cm^^km s"' 
toward W51 and W49N respectively. This is equivalent to as- 
suming a uniform HI optical depth in the gas components where 
we observe molecular absorption. Such an approximation un- 
derestimates the HI column density per unit velocity by no more 
than a factor 2. The total hydrogen column densities estimated 
with this method and given in Table A.l for W49N and W51 
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Table 4. HI, H2 and total hydrogen column densities. 









A'(H)" 


A'(H2)" 


* r, 

K 


Source 


(km s ) 


(km s ) 


(cm -) 


/ -''\ 

(cm -) 


(cm ^) 








xlO'^' 


xlO'^' 


xlO^' 




6.2 


5.6 


1.27 






W51 


11.8 


6.0 


0.58 


0.47 


2.3 








N{11)" 


A'(H2)^ 


> 


Source 


(km s-') 


(km s-') 


(cm-2) 


(cm-2) 


(cm-2) 








xlO^' 


xlO^' 


xlO-' 




J L.J 


-10 










-10.0 


4.0 


> 6.68 






4.0 


20.0 


5.87 






20.0 


29.3 


> 3.86 








-1.7 


11.3 


4.80 






G08.67-0.36 


11.3 


25.7 


5.32 




14 




25.7 


45.3 


> 9.62 








n.3 


25.0 


3.53 






GlO.62-0.38 


25.0 


30.6 


2.61 


15.8 


17 




30.6 


47.3 


5.54 








-2.5 


7.1 


0.54 






G34.3+0.1 


7.1 


20.0 


2.01 
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20.0 


34.1 


4.42 






W49N 


30.0 


50.0 


6.95 


1.6 


23 


50.0 


78.2 


7.23 


4.0 



" Gaussian decomposition from Koo (1997). The column densities 
have been scaled to a spin temperature of 100 K. 

* Inferred from the absorption profiles observed by Fish et al. (2003) 
with the VLA interferometer. A spin temperature of 100 K is as- 
sumed. 

Estimation from the CH emission lines observed by Rydbeck et al. 
(1976) toward W51 and W49N, the CH absorption lines observed 
by Gerin et al. (2010 in prep.) toward GlO.62-0.38, and the corre- 
lation between CH and H2, A'(CH)/A'(H2) = 4.3 x lO"** (Liszt & 
Lucas 2002). 

From models of the extinction at 2 jim. by Marshall et al. (2006). 



are smaller by only 50% than those inferred from IR extinction. 
We therefore estimate that the total H column density per ve- 
locity component on these two lines of sight does not exceed 
1.5 magnitude (or about 2.5 x 10^' cm"^). This result suggests 
that the HCO^ abundance in the different velocity components, 
defined as X(HCO+) = A^(HCO+)/A^h, ranges between 5x10""' 
and 5x10 toward both W49N and W51. For HCN, the scatter 
of abundances is also an order of magnitude among the com- 
ponents, with values ranging between 2x10"^ and 3x10"** for 
W49N. A similar scatter is also found for HNC with abundances 
~ 10 times smaller. 

In summary, since the column densities of the Gaussian com- 
ponents span less than two orders of magnitude while their 
linewidths span only a factor ~ 10 (between 0.3 and 3.4 km 
s"'), actual fluctuations of molecular abundances are therefore 
observed among the components. 

4.3. Physical properties of the gas seen in absorption 

The above average description is crude and is only meant to as- 
cribe an average column density to a velocity interval, ignor- 
ing velocity crowdings: the lines of sight toward distant star- 
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forming regions sample gas components with a broad distribu- 
tion of densities, velocity dispersions, and column densities. We 
have shown that the upper limit on the gas density is hh < 5 x lO"* 
cm"^ (cf. Sect. 13. lb and that the total column densities per veloc- 
ity components A^h are at most of the order of a few magnitudes, 
similar to those obtained along the high latitude lines of sight ob- 
served by Liszt & Lucas (2001). According to the definitions of 
Snow & McCall (2006), the gas sampled by these lines of sight 
is a mixture of diffuse (rtu < 500 cm"^ with a shielding from the 
UV field Ay < 1) and translucent gas ( 500 cm"^ < «h < 5000 
cm""* and a shielding 1 < Ay < 2). 

Allowing the gas density to range between 30 cm"^ and 
5 X 10^* cm""*, the total Nh per velocity component of a few 
magnitudes translates into sizescales ranging between ~ 0.2 and 
~ 30 pc. It is interesting that the coiTesponding range of velocity 
dispersions inferred from the linewidth-scale relation in the dif- 
fuse molecular gas (Falgarone et al. 2009) is ~ 0.2 - 5 km s"', 
very close to the observed range of component velocity disper- 
sions. It suggests that the mixture of diff'use and transluscent gas 
sampled by these lines of sight is entrained in the turbulent inter- 
stellar cascade. Although one velocity component is not neces- 
sarily a cloud, the linewidth-scale relation ascribes a sizescale to 
each velocity dispersion. Am ~ 1(^/1 pc)'^^ km s"'. Therefore, 
different molecular species present in a given velocity compo- 
nent of width Am are colocated within ~ 0.1 pc (if Am = 0.3 km 
s"') and 9 pc (if Am = 3 km s"'). It is the correlation of molec- 
ular species, observed at these grains, that we will now confront 
to the predictions of different numerical codes. 

In the following sections, we therefore study the chemistry 
of diffuse and translucent gas, and we compare the observed col- 
umn densities to the results of PDR (PhotoDissociation Region) 
and TDR (Turbulent Dissipation Region) models. 



and (3) the relaxation stages where the gas previously heated 
cools down to its original state. Because of the strong experi- 
mental (Mouri et al. 2007, Mouri & Hori 2009) and observa- 
tional (Crovisier 1981 ; Joncas et al. 1992 ; Miville-Deschenes 
et al. 2003, Haud & Kalberla 2007) constraints on uemm, only 
two main parameters govern the TDR model: the turbulent rate 
of strain that describes the coupling between the large and small 
scales of turbulence and the gas density. However, both parame- 
ters are not independent because of the constraint imposed by 
the observed turbulent energy available in the CNM, namely 
Sohs ~ 2 X IQ-^^ erg cm^^s"' (Falgarone 1998, Hily-Blant et 
al. 2008). As a result, the influence of the parameters on the 
TDR model, fully discussed in Godard et al. (2009), is complex. 

• The gas density riu affects the chemistry itself and has sev- 
eral impacts due to the dynamics and the geometry of the line 
of sight. For the range of parameters explored here, the dis- 
sipation of turbulent energy over one vortex is mainly due to 
ion-neutral friction and varies roughly as ; the vortex ra- 

— 1/2 

dius varies as ' , and their number along the line of sight 
as n^. 

• While the dissipation of turbulent energy over one vortex and 
the number of vortices along the line of sight are nearly in- 
dependent of the turbulent rate of strain a, their size and their 
lifetime vary roughly as a^^^^ and a respectively. 

For both models and as shown in Table IC.fl of Appendix ICl 
the nitrogen gas phase chemistry has been updated according to 
the recent measurements and rate review. However, neither the 
PDR nor the TDR chemical networks take into account the nitro- 
gen grain surface chemistry: this assumption will be discussed in 
Sect.lO 



5. Comparison of observations to model 
predictions 

5.1. Introduction to the PDR and TDR models 

The comparison of the PDR and TDR models has been discussed 
in detail in Godard et al. (2009). We recall here the main differ- 
ences between the two approaches. 

The PDR model is a 1 -dimensional chemical model in which 
a static slab of gas of uniform density and given thickness (or 
total column density of gas, noted A^h in the following) is illumi- 
nated by the ambient interstellar radiation field either on one side 
or on both sides (Le Petit et al. 2006). The computed column 
density of a molecular species to be compared to the observed 
values is therefore an integral performed over the slab thickness 
along a direction perpendicular to its surface. 

The TDR code is a 1 -dimensional model in which the chem- 
ical and thermal evolution of a turbulent dissipative burst - 
namely a magnetized vortex - is computed. The dynamics and 
lifetime of the magnetized vortices are controlled by the turbu- 
lent rate of strain a (in s"'), the gas density nu (in cm""*), and the 
maximal orthoradial velocity M^max (in cm s '). A random line 
of sight crossing several dissipative regions is then modelled by 
taking into account: the averaged turbulent dissipative rate ob- 
served in the CNM, and the long-lasting thermal and chemical 
relaxation stage that follows any dissipative burst. The resulting 
line of sight therefore intercepts three kinds of diffuse gas: (1) 
mainly the ambient medium (the coiTesponding filling factor is 
larger than 90 %) in which the chemistry is driven by the UV ra- 
diation field, (2) the active vortices where the gas is heated and 
the chemistry is enhanced by the dissipation of turbulent energy. 



5.2. Gas phase nitrogen chemical networks in the PDR and 
TDR models 

The networks of dominant reactions in the nitrogen gas phase 
chemistry are shown in Figs. lC.l l& IC.2l (AppendixlCli. While in 
both models, photodissociation largely dominates the destruc- 
tion mechanisms of CN-bearing molecules, there are several 
pathways that lead to their production. 

In the PDR model (Federman et al. 1994, Roger & Sternberg 
2005), two main formation routes are at work. The first involves 
the hydrogenation chain of carbon and the formation of CH and 
CH2, which in turn, leads to CN and HCN via the neutral-neutral 
reactions: 

CH H- N ^ CN H- H (5) 

and 



CH2 + N ^ HCN + H. 



(6) 



The second one involves the hydrogenation chain of nitrogen 
and the formation of NH and NH2, followed by the ion-neutral 
reaction chains 



NH^CN^^ ™; -^HCNH^ 



(7) 



and 



c+ 



NH2 HCN+ — * HCNH+, (8) 

the dissociative recombination of HCNH"^, and the photodisso- 
ciation of HCN andHNC. 

In the TDR model, because of both the differential rotation at 
the edge of the vortex and the ion-neutral drift, the gas is locally 
heated, and several key endoenergetic reactions such as 

C+ H- H2 ^ CH+ + H - AE/k = 4 640K (9) 
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open up. A warm nitrogen chemistry is triggered locally (Fig. 
IC.2I Appendix ICli. Although the temperatures reached in a vor- 
tex (~ 10^ K) are not suflicient to activate the reaction 



N H- H2 ^ NH -H H 



■AE/k = 16 606K 



(10) 



(whose rate becomes large enough only at 7000 K, Crawford 
& Williams 1997), the pathways leading to the production of 
N bearing species are deeply modified and their abundances 
increase considerably. The opening of the endoenergetic route 
C""-!- H2, that forms CH+ and then CH2+ and CH3+ by fast reac- 
tion with H2, triggers in turn the enhancement of the production 
of CN-bearing molecules through the ion-neutral reaction chains 



CH+ 



HCN^ 
HNC+ 



H2 



HCNH+ 



CH^ 



HCN+ 



HCNH+ 



(11) 



(12) 



and 



HCNH+ 



(13) 

followed by the dissociative recombination of HCNH^ and the 
photodissociation of HCN and HNC. Therefore, the key species 
here is CH^ . Its reaction with atomic nitrogen is able to signifi- 
cantly enhance the production of HCNH^ compared to that in the 
PDR models, leading to HCN, HNC, and CN. Similarly, CH3+ 
is responsible for the enhancement of HCO^ by its reaction with 
atomic oxygen (see the discussion in Godard et al. 2009). 

5.3. Comparison of PDR model predictions witli the 
observations 

To explore the role of the parameters in the Meudon PDR code 
(Le Petit et al. 2006) we computed several two-sided illuminated 
PDR models, as follows: 



• based on the results of Sect. 14.21 the total column density of 
gas in the slab is either A^h = 1-8 x 10^' or 3.6 x 10^' cm"^ 
(corresponding to 1 and 2 magnitudes respectively); 



• Two values of x are considered, x - ^ and = 3, in order 
to bracket the ambient UV radiation field found to vary be- 
tween the Solar Neighbourhood (~ 8.5 kpc from the Galactic 
centre) and the molecular ring (~ 4 kpc from the Galactic 
centre) by a factor ~ 3 (Moskalenko et al. 2006); 

• the gas density is set to «h = 10^, lO-', and 10"* cm"-'. 

Each model was computed assuming a cosmic ray ionisation rate 
^ = 3 X 10"'^ s"' (Dalgarno 2006). The results on the column 
densities and the column density ratios of HCO^, HCN, HNC, 
CN, and C2H are given in Fig. |6] which also displays, for com- 
parison, the ranges of observed values. 

This figure shows that the constraints provided by the CN- 
bearing molecules and by HCO^ cannot be reconciled in the 
framework of this model. On the one hand, both the range of 
column densities and column density ratios of CN, HCN, HNC, 
and C2H can be explained by PDR models for A^h > 1-8 x 10^^ 
cm"^ and a gas density larger than 300 cm"^. But on the other 
hand, for the same range of density, not only the predicted cor- 
relation between the column densities of CN and HNC is clearly 
non-linear and mismatch the observed value by a factor larger 
than 10, but the column densities of HCO^ are also underesti- 
mated by at least one order of magnitude. 

As expected, if the mean interstellar UV radiation field in- 
creases, the column densities of CN, HCN, HNC, and C2H de- 
crease because the photodissociation is the main destruction pro- 
cess of all these molecules. Oppositely, when x is multiplied by 
3, the temperature at the edge of the PDR is twice larger and 
the O -I- H^ charge exchange reaction (with an endothermicity 
AEjli of 227 K) is enhanced. Since this reaction is at the root of 
the formation of HCO^ in low density UV-dominated gas phase 
chemistry (Godard et al. 2009) and since the main destruction 
mechanism of HCO^ is its dissociative recombination (whose 
rate is independent ofx)^ the column density of HCO^ increases 
(see Fig.|6]l. All those behaviours holds for the ranges of density 
and radiative conditions explored in Fig.|6l 
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The dependence of the results on the cosmic ray ionisation 
rate ^ is slightly more complex. When ( is multiplied by 10, 
the abundances of and rise by one order of magnitude^. 
Thus, both the hydrogenation chains of nitrogen and oxygen are 
stimulated by the reactions 

N + H+^NH++H (14) 

and 

+ H+^0++H (15) 

respectively, while the hydrogenation chain of carbon, initiated 
by the slow radiative association 

C++H2^CH^+H (16) 

is independent of ^. Therefore when f is multiplied by 10: 

- the column density of HCO^ increases by a factor of 10 
whatever the gas density ; 

- since CN always originates from the carbon hydrogenation 
chain, its column density is independent of ( ; 

- since HCN and HNC originate from the carbon hydrogena- 
tion chain at large density and from the nitrogen hydrogena- 
tion chain at low density, their column densities are indepen- 
dent of ^ if hh > 300 cm""*, and increase by a factor of 2 if 
rin < 300 cm^^ 

Nevertheless, even at low density, the impact of the cosmic ray 
ionization rate is not sufficient enough to account for the range 
of column densities observed in Fig|6l 

5.4. Comparison of observations witli TDR model predictions 

As suggested by Godard et al. (2009) and by the present obser- 
vations, we explore the TDR models in the range of parameters: 
10"'^ < a < 10"'° s"' and 30 < nn < 500 cm"^ for an ambient 
radiation field characterized hy x = 1 • Since the influence of this 
parameter on the results of the PDR model (see Fig.|6]) seems to 
be small compared to the role of A^h and /ih, we did not explore 
a broader range of values foT x- The shielding from the ISRF is 
assumed uniform. Ay = 0.5 and the amount of gas sampled is 
fixed to A^H = 1.8 X 10^' cm . Note that, in the TDR model. Ay 
is the actual UV-shielding of the gas and is no longer identical to 
the total hydrogen column density sampled by the line of sight 
(as it is the case in the PDR model). 

Fig.|7]displays the comparison between the observed molec- 
ular column densities and the predictions of the TDR models. 
Since the results of the TDR model are simply proportional to 
A^H, the computed values for higher (or lower) values of A^h are 
easily inferred from the display. TDR models, without fine tun- 
ing of the parameters, are consistent with the data. Over a broad 
range of turbulent rates of strain 10"'^ s~' < a < 10"'" s"' and 
for a gas density 30 cm"^ < nn < 200 cm"^, the absolute col- 
umn densities and column density ratios of HCO^, C2H, CN, 
HCN, and HNC are reproduced. The case «h = 500 cm"^ is 
ruled out because the computed HCO^ column densities are too 
small compared to the observed values, even if the total hydro- 
gen column density A^h or the averaged turbulent dissipative rate 
s are multiplied by a factor of 10. For the CN-bearing species, all 
deriving from HCNH^ during the dissipation stage, the agree- 
ment with the observations holds over a broader range of gas 



^ and are formed via the reaction chains H2 — > H, — ^ H3 

CRP H2 

and He — > He^ — > respectively (CRP: cosmic ray particle). 



density : 30 cm"^ < nn < 500 cm"^. A trend, that was already 
found in Godard et al. (2009), is that the observed column den- 
sities are in better agreement with models of small turbulent rate 
of strain. 

6. Discussion 

6.1. Grain surface cliemistry 

As shown in the previous section, if only the gas phase chemistry 
is taken into account, the predictions of UV-dominated chemical 
models computed at moderate densities (riu < 300 cm"-') are un- 
able to explain the observed column densities of the CN-bearing 
species. This result is in agreement with the previous studies 
of Wagenblast et al. (1993), Crawford & Wilhams (1997), and 
Liszt & Lucas (2001) which show that the predictions of purely 
gas phase PDR models fail by a factor of 7 to account for the 
CN abundances observed in the low density material toward ( 
Per, and by factors of 30 and 40 to reproduce the NH abundances 
observed toward f Per and ( Oph respectively. 

In the present paper we explore the role of turbulent dissi- 
pation on the cyanide chemistry that operates via the hydrogena- 
tion chain of ionized carbon followed by the reaction chains (fTTT i. 
(fT2l i. and (fT3] l. An ongoing related work is the comparison of the 
TDR model predictions to the Herschel-HIFI observations of ni- 
trogen hydrides along W31C and W49N (Persson et al. 2010). 

But there is an alternative possibility to the pure gas phase 
chemistry. It has been proposed that the grain surface chemistry 
dominate all the other formation processes of the ammonia NH3, 
and is therefore responsibly for the production of NH and CN 
in the diffuse ISM (van Dishoeck & Black 1988). Assuming a 
NH3 grain formation rate coefficient of 

^G-5.4xl0-'«(-|-i^r cmV, (17) 
\miN) [amu]/ 

corresponding to a surface reaction efficiency twice lower than 
that for the H2 formation on grains, Wagenblast & Williams 
(1996) found that a UV dominated chemical model reproduces 
the CN and NH abundances observed toward ^ Per and f Oph. 
However, this approach is based on the assumption that if a ni- 
trogen atom strikes and sticks to a grain, ammonia is quickly 
formed and comes off the grain. The resulting rate coeficient, 
namely the cross section of the grain multiplied by the average 
speed of nitrogen atoms, clearly sets an upper limit to the effi- 
ciency of the grain surface chemistry. 

6.2. Tlie UV radiation field 

Although, as said above, the comparison of model predictions 
to observations is a difficult task because of the lack of informa- 
tion on the gas topology (density structure, division of the matter 
therefore shielding from the ambient UV field, homogeneity of 
the gas conditions among the different fragments) and the poorly 
known physical characteristics of what is called a "Gaussian- 
velocity component", we have three main constraints to guide 
the models. 

In the Gaussian-velocity components, (i) the gas density is 
lower than rtu = 5 x lO"* cm"^, where it is estimated, the 
column density of H2 is between 0.3 and 1.3 that of atomic hy- 
drogen (Table 4), and ( the average hydrogen column density 

* In a UV-dominated chemistry, the main destruction route of NH3 
is its photodissociation. The product NH2 then leads to the formation of 
NH and CN (see Fig.lCTTl. 
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1e+10 



1e+12 



1e+11 



1e+10 



1e+12 



N(HNC) (cm" 



N(HNC) (cm" 



Figure?. Observations compared to the predictions of TDR models. The data (open circles) are from Lucas & Liszt (2000), Liszt 
& Lucas (2001), Gerin et al. (2010), and this work. The TDR models (filled symbols) are computed for several densities : 30 
(triangles), 50 (squares), 100 (crosses), 200 (circles), and 500 (double crosses) cm"^. All models are computed for A^h = 1.8 x 10^' 
cm"^: note that if the total column density was larger, the model predictions would be simply multiplied by the same factor The 
TDR models are computed for Ay - 0.5 and a varying along each curve between 10 (top right) and 10 s ' (bottom left). 



per unit velocity interval ranges between 2 and 5 10^" cm"^/km 
s This implies that the total hydrogen column density cannot 
be much larger than 4 x 10^' cm"^ and that the shielding from 
the ambient radiation field is efficient (otherwise the H2 fraction 
would be small). According to PDR models, this imposes a nar- 
row combination of the total column density A^h and ambient 
UV field intensity x- We are therefore confident that the gross 
properties of the UV illumination of the gas in the "Gaussian 
velocity-components" are well estimated. 

6.3. The velocity field 

The fact that the widths of the HCO^ components are system- 
atically larger than those of the CN-bearing species suggests an 
involvement of the velocity field in the production and the evo- 
lution of these species. Observations in the visible domain by 
Lambert et al. (1990), Crane et al. (1995), and Pan et al. (2004, 
2005) already showed that the CH^ line profiles are broader and 
less Gaussian than those of CH, themselves broader than those of 
CN. It has been proposed that such behaviours are due to spatial 
confinement of the molecules whose production occurs hierar- 
chically into denser and colder environments (Crawford 1995). 



But the differences observed among the linewidths may also bear 
the signatures of the dynamical processes at work in the forma- 
tion of all those molecules convolved with their most different 
relaxation times (Godard et al. 2009). Last, the similarity be- 
tween the results obtained along Galactic lines of sight and in 
the Solar Neighbourhood also suggests a formation mechanism 
weakly dependent on the ambient UV field, i.e. possibly driven 
for instance by the dissipation of the ubiquitous turbulence. 

7. Summary & conclusions 

We have presented the analysis of single dish observations of 
mm absorption lines of HCO^, HNC, HCN, and CN toward re- 
mote star-forming regions in the Galactic plane. The density of 
the gas responsible for the absorption lines is constrained by the 
excitation temperature of the first rotational levels of HCO^ and, 
in two cases, by line emission of CO (and isotopes). An upper 
limit hh < 5 X lO^* cm""* is obtained. The estimation of the total 
column densities along the lines of sight, combined with their H2 
fraction, where known, gives an upper limit on the UV-shielding 
per velocity component. Ay < 2. These conditions suggest that 
the gas sampled along the lines of sight belongs to the diffuse 
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and/or translucent media. Because the decomposition of the ve- 
locity structure into Gaussian components has been carried out 
recursively for all the molecules and hyperfine components, in- 
tercomparisons of molecular abundances are feasible. The in- 
ferred velocity structure will be helpful to unravel the analysis 
of future observations with Herschel/HIFI. 

Our main result is the similarity between the column den- 
sities per velocity component derived in the Galactic plane gas 
(this work) and those derived in the Solar Neighbourhood in pre- 
vious works. It is unexpected because the spatial scales of the 
Galactic environment sampled in this work are much larger. This 
suggests that the physical and chemical processes involved in the 
formation and destruction of the analysed species are similar in 
diffuse or translucent environments over the whole Galactic disk. 
We also find that the HCO^ linewidths are broader than those of 
the CN-bearing molecules. 

We compare the observed column densities to two kinds of 
gas phase chemical models, namely the photodissociation region 
(PDR) and turbulent disspation regions (TDR) models and find 
that: (1) the observed column densities of HCO^ cannot be re- 
produced by PDR models of diffuse or translucent gas, (2) the 
column densities of CN-bearing species can be explained by 
PDR models applied to dense (hh > 300 cm and well shielded 
(Ay > 1) molecular gas, (3) the column densities and column 
density ratios of HCO+, CN, HCN, HNC, and C2H are repro- 
duced simultaneously by TDR models of diffuse and translucent 
gas, over a broad range of turbulent rates of strain (10"'^ s"' 
< a < 10"'° s"'), and gas densities (30 cm"^ < riu < 200 cm"^). 
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Appendix A: Gaussian decomposition and 
calculation of column densities 

Tabs. lA. l1 - IA.4l contains the results of the Gaussian decompostion procedure that 
we have applied to the spectra. The column densities, given in the last columns, 
are derived assuming a single excitation temperature T^x for all the levels of a 
given molecule as 

A' = e(rex)^--;^fl [rdv (A.l) 
gi, A„i 1 J J 

where vq, gu, gi and A„/ are the rest frequency, the upper and lower level de- 
generacies and the Einstein's coefficients of the observed transition, QiT^^) is 
the partition funtion, and c is the speed of light. We also remind that for a 
Gaussian profile of peak opacity tq and FWHM Av, the opacity integral is 
j Tdv = 1/2 V^/ In 2 To Av. For an excitation temperature of 2.73 K, Eq. lA.ll 
l)ecomes 

A'(HCO+) = 1.10 X 10'^ f'^'^'^ ^^-^^ 
N(HNC) = 1.78 X 10'^ f'^''^ ^^-^^ 
WCHCN) = 3.43X lO'M rdi; cm"^, (A.4) 
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and 



N{CN) = 6.96 X \0" J T dv cm"^ 



(A.5) 



using the HCO+ 7 = 0-1, HNC 7 = 0-1, HCN 7, F = 0, 1 - 1, 2, and CN 
J, Fl, F = 0, 1/2, 3/2 -1,1/2, 3/2 transitions respectively. 



Appendix B: Impact of the abscissa uncertainty on 
the multi-Gaussian decomposition procedure 

In a spectrum, the possible velocity substructures are systematically erased due 
to the finite velocity resolution Sv. This could be modelled as an uncertainty on 
the velocity position of each point. Unfortunately the errors on the abscissa are 
rarely included in non-Unear fitting procedures because the system is consider- 
ably heavier to solve and because it often prevents convergence. To evaluate the 
resulting uncertainties on the fit pai'ameters, namely the central opacity, the ve- 
locity centroid, and the FWHM, we apply the fitting procedure on 3000 synthetic 
spectra of FWHM Am varying between between 0.3 and 3.4 km s"' sampled with 
the finite spectral resolution of the observations 6v ~ 0.13 km s"'. A noise is 
added to the .v-coordinates of all the spectral points. The rms of all the measured 
linewidths Au' is found to scale as 



d{Av') 
Av 



0.45 



Iam^ 



(B.l) 



and decreases from 0.04 to 0.01 km s"' as the true linewidth increases from 0.3 to 
3.4 km s"' . These uncertainties are smaller than (or comparable to) those inferred 
from the fitting procedure and the resulting errors on the column densities are at 
most 12%. In comparison the resulting eiTors on central opacities and velocity 
centroids are negligible. 



Appendix C: Cyanides chemical network 

Figs. IC. II & IC.2I show the main production and destruction pathways of the 
hydrogenation chains of carbon, nitrogen, and cyano, resulting from the PDR 
(«H = 50 cm-3, Av = 0.4) and TDR (tin = 50 cm-\ Ay = 0.4, o = 10"" s"') 
models respectively. These figures ai'e simplified: for each species, only the re- 
actions which altogether contribute at least to 70 percent of the total destruction 
and formation rate are displayed. There is one major difference between these 
networks: in a UV-dominated chemical model, the cyanide chemistry is initiated 
by: 

CH + N^CN + H, (C.l) 



and 



CH, + N HCN + H, 



NH + C+ ^ CN+ + H; 



(C.2) 



(C.3) 



while in a chemistry driven by turbulent dissipation, the hydrogenation chain of 
cyano is triggered by the ion-neutral reactions: 



CH+ + N ^ CN+ + H, 
CH^ + N ^ HCN+ + H, 



and 



CH+ + N - 



HCNH+ + H. 



(C.4) 
(C.5) 

(C.6) 



Since the pathways displayed in Figs. lC.ll & IC.2l depend on the chemical rates, 
and since the nitrogen and cyanide chemistry are still poorly known, we list the 
chemical rates we have adopted in our models for several reactions in Table lC. 11 




-0.5 0.0 0.5 1.0 1.5 2.0 2.5 




-0.5 I I 

-0.5 0.0 0.5 1.0 1.5 2.0 

''^HCO'^(O-I) 

Figure 2. HCO+ (1-2) opacities as a function of HCO+ (0-1) 
opacities. The lines are the results of linear regressions whose 
only parameter is the excitation temperature Tex- 
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Table A.l. HCO^ (0-1) absorption line analysis results. 







Au 


To 


A'(HCO+) 




Source 


(Icm s-') 












cm"") 


(10^" cm-2) 




-4. Q9 


1.13 ± 


0.15 


0.09 ± 


0.01 


0.13 


± 0.03 




G08.67-0.36 


9.82 


0.95 + 


0.12 


1.60 ± 


0.28 


2.02 


± 0.58 






11.71 


1.40 + 


0.17 


1.65 ± 


0.28 


3.08 


± 0.90 






14.52 


2.00 ± 


0.18 


1.41 ± 


0.18 


3.74 


± 0.83 






17.76 


1.17 ± 


0.22 


0.87 ± 


0.20 


1.34 


± 0.56 






20.00 


2.55 ± 


0.63 


1.08 ± 


0.16 


3.65 


± 1.45 






22.00 


1.02 ± 


0.58 


0.24 ± 


0.17 


0.32 


± 0.41 




GlO.62-0.38 


11.62 


1.10± 


0.12 


0.43 ± 


0.04 


0.62 


± 0.12 






16.38 


2.35 ± 


0.10 


2.23 ± 


0.14 


6.84 


± 0.74 






19.00 


2.13 ± 


0.20 


1.50 ± 


0.09 


4.18 


±0.64 






22.20 






1 S7 + 
i . J / m 




S S7 
J.J / 


+ 1 nn 








1.15 ± 


0.45 


0.82 ± 


0.18 


1.23 


±0.75 






24.98 


1.65 ± 


0.10 


0.70 ± 


0.06 


1.51 


± 0.21 






27.50 


2.58 ± 


0.17 


1.61 ± 


0.11 


5.46 


± 0.74 






29.00 


1.53 ± 


0.52 


0.78 ± 


0.20 


1.56 


± 0.93 






30.50 


1.52 ± 


0.32 


2.09 ± 


0.27 


4.13 


± 1.39 






40.91 


2.60 ± 


0.10 


2.33 ± 


0.14 


7.93 


± 0.78 






44.94 


1.22 ± 


0.05 


1.30 ± 


0.07 


2.06 


± 0.21 






46.40 


fi fin + 

U.UU It 


n in 


n 9fi + 


n rizi 


n 9n 


+ n 07 




G34.3-I-0.1 


9.80 


0.53 ± 


0.22 


0.38 ± 


0.12 


0.27 


±0.19 






10.60 


0.57 ± 


0.12 


1.86 ± 


0.61 


1.38 


± 0.72 






11.70 


1.37 ± 


0.18 


1.82 ± 


0.39 


3.34 


± 1.14 






14.07 


0.37 ± 


0.08 


0.77 ± 


0.23 


0.38 


± 0.20 






27.00 


1.02 ± 


0.13 


1.19± 


0.20 


1.61 


± 0.47 






27.75 


1.83 ± 


0.32 


0.45 ± 


0.09 


1.10 


± 0.41 






48.30 


2.23 ± 


0.20 


0.94 ± 


0.10 


2.80 


±0.54 




W49N 


33.20 


0.65 ± 


0.03 


0.44 ± 


0.02 


0.34 


± 0.02 


3.0 




33.70 


2.23 ± 


0.03 


0.63 ± 


0.01 


1.67 


± 0.02 


10.5 




36.50 


1.67 ± 


0.12 


0.29 ± 


0.01 


0.60 


± 0.03 


7.8 




38.50 


1.75 ± 


0.08 


0.91 ± 


0.01 


1.88 


±0.03 


8.2 




39.70 


1.47 ± 


0.02 


2.01 ± 


0.03 


3.47 


± 0.05 


6.9 




50.92 


1.65 ± 


0.10 


0.35 ± 


0.01 


0.68 


± 0.02 


7.8 




53.40 


1 7^ + 


0.10 


n ^n + 


0.02 


0.43 


+ 0.03 


5.8 




"^d on 


3.38 ± 


0.15 


0.51 ± 


0.01 


2.03 


± 0.04 






56.70 


2.35 ± 


0.07 


0.69 ± 


0.01 


1.91 


± 0.03 


11.0 




59.55 


2.35 ± 


0.05 


1.89 ± 


0.02 


5.24 


± 0.07 


11.0 




60.80 


0.88 ± 


0.05 


0.77 ± 


0.05 


0.80 


± 0.06 


4.1 




62.65 


2.95 ± 


0.03 


1.48 ± 


0.01 


5.16 


± 0.05 


13.9 




66.21 


1.85 ± 


0.07 


0.19 ± 


0.01 


0.41 


± 0.01 


8.7 




68.58 


1.80 ± 


0.10 


0.24 ± 


0.01 


0.52 


± 0.01 


8.5 




71.11 


2.25 ± 


0.20 


0.07 ± 


0.00 


0.19 


± 0.01 


10.6 


W51 


3.80 


0.47 ± 


0.12 


0.10 ± 


0.02 


0.06 


± 0.03 


1.0 




5.10 


1.23 ± 


0.05 


0.72 ± 


0.02 


1.14 


± 0.08 


2.7 




6.25 


0.65 ± 


0.03 


1.05 ± 


0.06 


0.87 


± 0.10 


1.4 




6.85 


2.01 ± 


0.05 


0.70 ± 


0.02 


1.80 


±0.10 


4.4 




12.71 


1.67 ± 


0.18 


0.10 ± 


0.01 


0.20 


±0.04 


3.7 


b 




^max 


j Tdv 


A'(HCO+) 




Source 


(km s-') 


(km 1 


r') 


(1cm i 


r') 


(10'2 


cm -) 




GlO.62-0.38 


13.6 


15. 


2 


0.77 ± 


0.08 


0.94 


±0.10 






31.0 


33.0 


3.77 ± 


0.31 


4.63 


±0.38 






33.3 


34.7 


1.81 ± 


0.15 


2.23 


± 0.18 






34.7 


36.0 


1.74 ± 


0.16 


2.13 


±0.19 






36.0 


37.5 


3.10 + 


0.29 


3.82 


±0.36 






37.5 


40.0 


> 6.27 


> 


7.71 




W49N 


41.0 


47.0 


0.90 ± 


0.03 


1.00 


±0.04 





" Multi-Gaussian decomposition. The column densities are derived assuming the excitation temperatures given in Table|3] 
'' Failure of the multi-gaussian decomposition. Column densities are inferred from the integral of the opacity assuming the excitation tempera- 
tures given in Table [3] 

Determined from the average total hydrogen column density per unit velocity A'h/Av = 2.2 and 4.7x10^" cm"^/km s"' toward W51 and 
W49N respectively. 



B. Godard et al: Molecular absorption lines toward star-forming regions 



15 




Cavbon Carbon & Nitrogen Nitrogen 

Figure C.l. Chemical network of a UV-dominated chemistry: nu - 50 cm"^ and Ai/ - 0.4. This figure is simplified: for each species, 
only the reactions which altogether contribute at least to 70 percent of the total destruction and formation rate are displayed. The 
red arrow show the endoenergetic reactions with the energy involved. 




Carbon Carbon & Nitrogen Nitrogen 

Figure C.2. Same as Fig. IC.ll for a turbulence-dominated chemistry: rtu = 50 cm"^, Ay = 0.4 and a = 10 " s"'. For the sake of 
simplicity, the main destruction route of N (N H- CHj — > HCN^ + H) and the main formation pathway of (photodissociation of 
NO^) are not displayed here: those two species are therefore highlighted. 
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Table A.2. HNC (0-1) absorption line analysis products. 





"0 


Av 


To 


A'(HNC) 


Source 


(km s-') 


(km ! 


r') 






(10'2 


cm -) 


GlO.62-0.38 


11.62 


1.72 ± 


0.22 


0.19 + 0.02 


0.72 


± 0.17 




14.42 


1.23 ± 


0.42 


0.09 H 


t0.02 


0.24 


± 0.14 




16.38 


1.30 ± 


0.05 


1.43 H 


t0.07 


4.02 


±0.33 




19.16 


2.01 ± 


0.15 


0.50 H 


t0.03 


2.19 


± 0.27 




22.20 


2.15 ± 


0.28 


0.38 H 


t0.06 


1.76 


±0.51 




22.95 


1.17 ± 


0.13 


1.21 H 


t0.08 


3.06 


±0.57 




24.98 


2.05 ± 


0.35 


U. i J □ 


t U.UZ 


0.60 


±0.18 




27.80 


1.50 ± 


0.07 


1.18 H 


t0.08 


3.82 


± 0.44 




29.10 




U. 1 o 


1.16H 


t 0.10 




-t- n 






1.13 ± 


0.07 


1.23 H 


t 0.10 


3.02 


± 0.41 




31.90 


2.31 ± 


0.37 


n ^ 
u.^^ □ 


t U.Uj 


2.28 


± 0.50 




34.00 


1.08 ± 


0.22 


0.32 H 


t0.04 


0.76 


± 0.25 




35.35 


1 09 -t- 


U. 1 J 


0.53 H 


t0.05 


1 1 fx 

i . ID 


± U.ZD 






1.22 ± 


0.10 


0.81 H 


t0.05 


2.12 


±0.31 




38.80 


2.16 + 


0.18 


0.72 H 


t0.03 


3.38 


±0.45 




40.91 


1.30 ± 


0.05 


2.16 H 


t 0.14 


6.10 


±0.63 




44.94 


1.30 ± 


0.10 


0.34 H 


t0.03 


0.95 


± 0.15 


G34.3-I-0.1 


10.60 


0.58 + 


0.13 


0.15 H 


t 0.03 


0.19 


± 0.08 




11.80 


1.45 ± 


0.15 


0.28 H 


t0.02 


0.85 


± 0.15 




27.00 


1.05 ± 


0.07 


0.48 H 


t0.03 


1.06 


± 0.13 


W49N 


33.50 


1.77 ± 


0.17 


0.12 H 


tO.Ol 


0.39 


± 0.06 




38.50 


1.68 ± 


0.10 


0.22 H 


tO.Ol 


0.70 


± 0.08 




39.70 


0.95 ± 


0.02 


1.46 H 


t0.04 


2.63 


± 0.11 




50.92 


0.53 ± 


0.15 


0.06 H 


tO.Ol 


0.06 


±0.03 




54.00 


2.85 ± 


0.32 


0.09 H 


tO.Ol 


0.47 


± 0.09 




56.70 


1.78 ± 


0.18 


0.10 H 


t 0.01 


0.32 


± 0.06 




59.40 


2.01 ± 


0.08 


0.44 H 


t 0.01 


1.67 


± 0.11 




60.90 


0.97 ± 


0.13 


0.08 H 


tO.Ol 


0.15 


± 0.04 




62.60 


2.15 ± 


0.08 


0.39 H 


tO.Ol 


1.58 


± 0.11 


b 




^max 


Irdv 


A'(HNC) 


Source 


(kms-') 


(Icm ! 


r') 


(km 


s-') 


(10'2 


cm -) 


GlO.62-0.38 


46.0 


47.0 


<0.03 


< 1 


0.07 


G34.3-H0.1 


9.0 


10. 


1 


< 0.03 


< 1 


0.06 




LJ.J 


14.7 


< 0.04 


< 1 


0.07 




46.0 


50.0 


<0.06 


< 1 


0.11 


W49N 


35.0 


37.0 


0.09 = 


t0.02 


0.16 


±0.04 




41.0 


47.0 


0.11 = 


t0.04 


0.19 


±0.06 




65.0 


67.3 


0.10 H 


t0.02 


0.18 


± 0.04 




67.3 


69.7 


0.07 H 


t0.02 


0.12 


± 0.04 




69.7 


72.5 


<0.02 


< i 


0.04 


W51 


3.0 


5.5 


0.09 H 


t 0.05 


0.17 


± 0.10 




5.5 


9.5 


0.09 H 


t 0.06 


0.18 


± 0.12 




11.0 


14.0 


<0.06 


< i 


0.11 



" Same as TablelAm 
* Same as TablelAll 
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Table A.3. HCN (0-1) absorption line analysis products. 







Av 




To 


A'(HCN) 




(km s"') 


(km ; 








(10'2 


cm -) 


GlO.62-0.38 


44.97 


0.92 ± 


0.03 


1.45 


±0.07 


5.40 


±0.45 


G34.3-H0.1 


10.10 


1.00 + 


0.15 


0.16 


± 0.02 


0.66 


±0.18 




10.60 


0.32 ± 


0.03 


0.40 


± 0.04 


0.53 


± 0.12 




11.80 


1.22 ± 


0.05 


0.85 


±0.03 


4.27 


±0.33 




27.00 


0.57 ± 


0.02 


1.47 


±0.06 


3.44 


± 0.24 




27.75 


0.55 ± 


0.05 


0.35 


±0.03 


0.81 


± 0.14 


W49N 


33.40 




U.UJ 


0.42 


± 0.01 


Z.ZU 


± U.UD 






1.22 ± 


0.03 


U.U J 


IC U.U i 


2.91 


± 0.06 




39.70 


0.87 ± 


0.02 


3.94 


±0.06 


12.51 


± 0.18 




51.00 


1.00 ± 


0.17 


0.11 


±0.01 


0.40 


± 0.11 




54.00 


2.50 ± 


0.20 


0.24 


±0.01 


2.23 


± 0.21 




56.60 


1.83 ± 


0.15 


0.28 


± 0.01 


1.85 


± 0.22 




59.40 


1.50 ± 


0.05 


1.35 


±0.02 


7.44 


±0.10 




60.90 


0.83 ± 


0.05 


0.64 


±0.02 


1.96 


±0.05 




62.55 


1.83 ± 


0.05 


1.15 


± 0.01 


7.73 


± 0.09 


W51 


5.00 


0.93 ± 


0.12 


0.17 


± 0.02 


0.63 


± 0.14 




6.80 


1.27 ± 


0.12 


0.20 


± 0.02 


1.02 


± 0.17 


b 






i 


rdu 


A'(HCN) 


Source 


(kms-') 


(km ! 


r') 


(km s-i) 


(10'2 


cm-2) 


W49N 


69.7 


72.5 


0.03 


±0.03 


0.18 


±0.15 


W51 


11.0 


14.0 


< 


0.04 


<0.27 



" Same as TablelATTI 
Same as TablelAll 
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Table A.4. CN (0-1) absorption line analysis products. 







Av 




To 


yv(CN) 


Source 


(kms-') 


(kms-') 






(10" 


cm--) 


GlO.62-0.38 


11.62 


0.92 ±0.18 


0.16 


±0.02 


1.19 


±0.35 




16.30 


1.20 ±0.03 


1.16 


± 0.03 


11.10 


1 ± 0.63 




19.16 


1.03 ±0.10 


0.47 


±0.02 


3.91 


± 0.59 




22.20 


1.83 ±0.00 


0.32 


±0.02 


4.75 


±0.35 




22.90 


0.87 ± 0.05 


1.16 


±0.05 


8.04 


±0.87 




27.80 


1.47 ±0.05 


1.23 


±0.04 


14.34 


■ ± 1.05 




29.10 


0.78 ± 0.07 


0.99 


±0.07 


A 90 


-1-0 09 










It u.uu 


12.28 


; ± 1.05 




31.80 


2.10 ±0.28 


0.38 


± 0.02 




-1-119 
It I. IZ, 






1 97 + n 78 

i ..iZ It \J.jLO 


U. i V 


+ fl 09 
It u.uz 


1.88 


±0.63 




35.35 


0.68 ± 0.07 


0.44 


± 0.03 


2.38 


± 0.38 




36.61 


1.40 ±0.13 


0.50 


± 0.02 


5.56 


± 0.77 




38.80 


1.63 ±0.10 


0.66 


±0.02 


8.65 


±0.84 




40.75 


1.35 ±0.03 


1.70 


±0.05 


18.30 


1 ± 0.95 




44.94 


0.83 ± 0.08 


0.22 


±0.04 


1.43 


±0.41 


G34.3-H0.1 


27.00 


0.68 ± 0.12 


0.86 


± 0.17 


1.79 


± 0.66 


W49N 


33.30 


1.83 ± 0.12 


0.13 


± 0.01 


1.82 


± 0.22 




38.40 


1.77 ± 0.08 


0.23 


± 0.01 


3.02 


± 0.27 




39.60 


0.85 ± 0.02 


1.18 


±0.03 


7.43 


±0.29 




59.20 


1.98 ±0.08 


0.27 


±0.01 


3.98 


±0.34 




62.40 


1.62 ± 0.07 


0.32 


±0.01 


3.77 


±0.31 


h 






.1 


rdv 


A'(CN) 


Source 


(kms-') 


(kms-') 


(kms-') 


(10" 


cm--) 


GlO.62-0.38 


13.6 


15.2 


0.18 


±0.06 


0.85 


±0.29 




24.0 


26.0 


0.26 


± 0.07 


1.29 


± 0.34 




46.0 


47.0 


< 


0.04 


< 


0.21 


G34.3-H0.1 


9.0 


10.1 


0.20 


±0.13 


1.54 


± 1.01 




10.1 


11.0 


< 


0.11 


< 


0.84 




11.0 


13.0 


0.48 


±0.18 


3.63 


± 0.14 




13.3 


14.7 


0.15 


±0.14 


1.12 


± 1.05 




4fi 
^u.u 




< 


0.19 


< 


1.46 


W49N 


35.0 


37.0 


< 


0.03 


< 


0.13 




41.0 


47.0 


< 


0.05 


< 


0.22 




50.0 


52.0 


< 


0.03 


< 


0.13 




52.5 


55.5 


0.09 


±0.04 


0.41 


± 0.17 




55.5 


57.9 


0.07 


±0.03 


0.32 


± 0.14 




60.3 


61.3 


0.09 


±0.02 


0.41 


± 0.10 


W51 


3.0 


5.5 


0.13 


±0.08 


0.97 


± 0.56 




5.5 


9.5 


< 


0.09 


< 


0.69 




11.0 


14.0 


< 


0.08 


< 


0.59 



" Same as TablelAll 
* Same as TablelAJl 
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Table C.l. Rates k of the main reactions of the cyanide chemistry. The photoreaction rates writte k - a ■ exp(-yAv) while all the 
other reaction rates are derived ask = a ■ (T /300Kfexp(-y/T) where T is the effective temperature (Zsargo & Federman 2003). 
Numbers in parenthesis are power of 10. 





Reaction 




a 


yS 


r 


ref 


CN 


+ hv 


^ C 


+ N 


4.64 (-10) 


0.00 


3.6 


a 


HCN 


+ hv 


^ CN 


+ H 


5.70 (-10) 


0.00 


2.6 


a 


HNC 


+ hv 


CN 


+ H 


5.70 (-10) 


0.00 


2.6 


a 


HCNH+ 


+ e" 


^ HCN 


+ H 


9.45 (-07) 


-0.65 


0.0 


b,c 


HCNH* 


+ e" 


^HNC 


+ H 


9.45 (-07) 


-0.65 


0.0 


b,c 


HCNH+ 


+ e" 


^CN 


+ 2H 


9.10 (-08) 


-0.65 


0.0 


b,c 


HCNH+ 


+ e" 


^CN 


+ H2 


1.80 (-08) 


-0.65 


0.0 


b,c 


H2NC+ 


+ e" 


^ HNC 


+ H 


1.75 (-07) 


-0.50 


0.0 


d 


H2NC+ 


+ e" 


^ CN 


+ H2 


1.75 (-07) 


-0.50 


0.0 


d 


HNC+ 


+ e" 


^ CN 


+ H 


3.90 (-07) 


-1.00 


0.0 


e,f 


HCN+ 


+ e" 


CN 


+ H 


3.90 (-07) 


-1.00 


0.0 


e 


CH3+ 


+ N 


^ HCN+ 


+ H2 


2.20 (-11) 


0.00 


0.0 


g 


CH3+ 


+ N 


^ HCNH+ 


+ H 


2.20 (-11) 


0.00 


0.0 


g 


CH3^ 


-l-N 


^ HzNC-^ 


+ H 


2.20 (-11) 


0.00 


0.0 


g 


CH2* 


-I-N 


^HCN-" 


+ H 


4.70 (-10) 


0.00 


0.0 




CH2* 


+ N 


^HNC* 


+ H 


4.70 (-10) 


0.00 


0.0 




CH+ 


+ N 


^ CN+ 


+ H 


1.90 (-10) 


0.00 


0.0 


h 


CH2 


+ N 


^ HCN 


+ H 


3.95 (-11) 


0.17 


0.0 


i 


CH 


+ N 


^ CN 


+ H 


2.00 (-10) 


0.00 


0.0 


j 


C+ 


+ NH 


^ CN+ 


+ H 


7.80 (-10) 


0.00 


0.0 


C+ 


+ NH2 


^ HCN+ 


+ H 


1.10 (-09) 


0.00 


0.0 




HCN+ 


+ H2 


^ HCNH+ 


+ H 


9.00 (-10) 


0.00 


0.0 


k 


HNC* 


+ H2 


^ HCNH* 


+ H 


9.00 (-10) 


0.00 


0.0 


k 



" Robergee^a/. (1991). 

* Florescu-Mitchell & Mitchell (2006). 
Semaniak aZ. (2001). 

Adams & Smith (1988). 
" Sheehane«aZ. (1999). 
^ Talbie/a/. (2000). 
« Fehsenf eld (1976). 

* Viggiano ef a/. (1980). 
' Smith et al. (2004). 

^ BrownswordeJaZ. (1996). 

* Hunti-ess (1977). 
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